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ABSTRACT 

The wide-band energy spectrum of NGC 6240 over the range 0.5-200 keV is investi- 
gated using the RXTE and ASCA data. The RXTE data provide the spectrum beyond 
the ASCA range (0.5-10 keV) with significant detection of signals up to 20 keV and 
the upper limits above 20 keV. The spectrum above 10 keV is found to be very flat. 
A strong iron-K emission line discovered in the previous ASCA observation is also 
confirmed with the RXTE PCA. These results provide further evidence for the dom- 
inance of a reflection component, i.e. emission from cool material illuminated by an 
AGN. By fitting the spectra obtained with RXTE and ASCA simultaneously, we sat- 
isfactorily modeled the AGN spectrum with a Compton reflection component and 
probably a transmitted AGN component penetrating through a thick absorber. The 
X-ray luminosity of the AGN is estimated to be in the range 4 x 10 43 — 6 x 10 44 
ergs/s in the range 2-10 keV, which categorizes NGC 6240 among the most luminous 
Seyfert nuclei. The ratio of the 2-10 keV X-ray luminosity to the infrared luminosity, 
Lx(2 — 10 keV)/-LiR, is 0.01 - 0.1, which implies a quite substantial, if not dominant, 
contribution of AGN to the infrared luminosity. 
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1 INTRODUCTION 

The IRAS survey (Neugebauer et al. 1984) discovered many 
ultraluminous infrared galaxies (ULIRGs) that emit the 
bulk of their energy in infrared (IR) photons. Since their 
bolometric luminosity and the number density are as high 
as those of quasars, ULIRGs are among the most energetic 
objects in the universe. The most fundamental problem yet 
to be solved is the energy source of the extremely intense 
infrared emission. 

NGC 6240, a gravitationally interacting system with a 
complex optical morphology (Fosbury & Wall 1979; Fried 
& Schulz 1983), is a very interesting example of ULIRG. Its 
bolometric luminosity reaches 2.4x 10 12 Lq (Weight, Joseph, 
& Meikle 1984; z=0.0245 and H = 50 km s^Mpc -1 are 
assumed). NGC 6240 is outstanding in several respects. Its 
H 2 1 -> 05(1) at 2.121/mi and [Fell] 1.644/im line luminosi- 
ties and the ratio of H2 to bolometric luminosities are the 
largest currently known (e.g., van der Werf et al. 1993). Fur- 
ther, its stellar velocity dispersion of 360 km/s is among the 
highest values ever found in a galaxy centre (e.g., Doyon et 
al. 1994). 



The energy source of the huge IR luminosity is con- 
troversial. Many IR spectroscopic studies (e.g. Genzel et 
al. 1998; Ridgway et al. 1994; Rieke et al. 1985; Weight et 
al. 1984) have suggested that main energy source of the IR 
emission is starburst activity, which is presumably a super- 
starburst induced by a merger of two galaxies (Joseph & 
Wright 1985; Chevalier & Clegg 1985). The ground-based 
optical spectrum can be classified as LINER, and is inter- 
preted as a result of shock heating (Heckman et al. 1987). 
On the other hand, a significant contribution from an ac- 
tive galactic nucleus (AGN) similar to Seyfert galaxies was 
also discovered from IR spectroscopy (DePoy et al. 1986). 
Another hint of an AGN in NGC 6240 is the presence of com- 
pact bright radio cores (Carral et al. 1990 but see Colbert et 
al. 1994). HST discovered a core that is excited higher than 
LINER (Rafanelli et al. 1997). 

X-ray observation provides an important tool for inves- 
tigating both the starburst and AGN activity. The ROSAT 
(Triimper 1990) observations showed that NGC 6240 is fairly 
bright in the soft X-ray band below 2 keV with a luminosity 
larger than 5 x 10 42 ergs s" 1 (0.1-2.0 keV), and the bulk of 
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the X-ray emission is extended in a scale of ~ 25" (Schulz 
et al. 1998; Komossa, Schulz, & Greiner 1998; Iwasawa & 
Comastri 1998). Detailed spectroscopic studies with ASCA 
(Tanaka, Inoue, & Holt 1994) data of NGC 6240 have shown 
that the soft X-ray spectrum can be explained with two 
thermal components, a cooler component with a tempera- 
ture of 0.2-0.6 keV and a hotter component of ~ 1 keV with 
an excess absorption of ~ 10 22 cm -2 (Iwasawa & Comastri 
1998). These results show that these soft X-rays are most 
likely originated from thermal processes, which may arise 
from starburst activities (e.g. Heckman et al. 1987). 

On the other hand, in the 3-10 keV band, another very 
hard continuum with a strong iron-K emission feature was 
observed with ASCA (Mitsuda 1995; Iwasawa & Comastri 
1998). The observed spectral feature can be accounted for 
in terms of Compton reflection from optically thick material 
(e.g. Lightman & White 1988; George, Nandra & Fabian 
1990), and is generally accepted as evidence for the pres- 
ence of an AGN in NGC 6240. The emission-line profile fur- 
ther indicates that a part of the reflector is highly ionized 
(Mitsuda 1995; Iwasawa & Comastri 1998), and both the 
ROSAT and ASCA spectra were modeled with a reflection 
component from warm material surrounding the AGN (Net- 
zer, Turner, & George 1998; Komossa et al. 1998). 

As described above, many observational facts of 
NGC 6240, in particular X-ray spectroscopy, supports the 
presence of an AGN in NGC 6240, which may account for a 
significant fraction of the huge IR luminosity. However, the 
intrinsic power of the AGN is still uncertain, The hidden 
AGN could be visible as a strongly absorbed X-ray contin- 
uum above 10 keV, penetrating through a thick layer of the 
obscuring matter. 

Very recently, Vignati et al. (1999) has published the 
BeppoSAX results of NGC 6240, concluding the detection 
of the direct X-rays from an AGN. As described in this 
paper, we performed an independent study of NGC 6240 
with RXTE which carries the PCA ( Jahoda et al. 1996) and 
HEXTE (Rothschild et al. 1998) covering 2-250 keV. 

In the spectral analysis, we also utilized the ASCA data 
covering 0.5-10 keV energy band with two Gas Imaging 
Spectrometers (GIS: Ohashi et al. 1996) and two X-ray CCD 
cameras (SIS: Burke et al. 1991). Simultaneous use of the 
ASCA data provides more constraints in modeling the spec- 
trum. Moreover, because RXTE has no imaging capability, 
the RXTE spectrum may be subject to contamination by 
nearby sources. In fact, we find from the ASCA image that 
it is the case, and the ASCA data are used to correct the 
RXTE spectrum. 

We have examined various spectral models, and con- 
clude the presence of a high-luminosity AGN in NGC 6240. 
Our results are in essential agreement with those of Vignati 
et al. (1999). 



2 RXTE DATA 

2.1 Observation and data reduction 

The RXTE observation of NGC 6240 was performed on 1997 
Nov 9-11. The data were reduced with Ftools 4.1 and 4.2, 
We discarded the PCA and HEXTE data that were taken 
when the centre of field of view was within 10° of the local 



horizon, and when it was off NGC 6240 by more than 0.01 
degree. We also selected the time period when all the five 
PCUs of PCA were on. The total on-source exposure time 
after the data screening process is 31.4 ksec and 10.3 ksec 
for the PCA and HEXTE, respectively. All the five PCUs' 
top layers of PCA are used for the analysis, while one of 8 
detectors of HEXTE (#3 in cluster B) is excluded because of 
its failure. We found no significant time variability from the 
PCA background-subtracted light curve in the entire energy 
band (3-20 keV) as well as in several different energy bands. 

2.2 Derivation of PCA and HEXTE spectra 

We estimated the PCA background using the current stan- 
dard method. ^ We attempted to estimate systematic error 
associated with the background subtraction by comparing 
the estimated background with the on-source data in the 
range above 30 keV where the source signal in PCA should 
be negligible due to low detection efficiency. The count rates 
of both agreed with each other within 0.01%. Therefore, we 
did not introduce any systematic error to the PCA back- 
ground. 

We estimated the HEXTE background using the data 
taken during the off-target pointings (Rothschild et al. 
1998). In order to estimate the HEXTE-background system- 
atics, we accumulated spectra from the +1.5° off-pointings 
and the -1.5° off-pointings, separately, and found that the 
difference between the two background spectra was within 
±2%. Therefore, we introduced a 2% systematic error to the 
HEXTE background spectrum. 

Figure |l|a illustrates the PCA and HEXTE energy spec- 
tra thus obtained. Significant source signals are detected up 
to ~25keV in the PCA data. However, the PCA data be- 
low 4 keV and above 20 keV were discarded for the spec- 
tral analysis below, due to the uncertainty of the energy 
response reported by Gierlinski et al. (1999) and Jahoda et 
al. (1996). The HEXTE spectrum shows significant signals 
up to 20 keV, while only upper limits are obtained above 
20 keV. 

2.3 Spectral analysis 

In order to see the overall spectral characteristics, we first 
fit the spectra with a simple power-law model modified by 
Galactic absorption, i.e. f(E) = A"£r r e~ <T p hJVH G - 1 , where 
E is X-ray energy, f(E) is photon flux given in units of pho- 
tons s _1 cm~ 2 keV" 1 , F is photon index, K is the nor- 
malization factor, <r p h is the cross-section of photoelectric 
absorption given by Morrison & McCammon (1983), and 
Ah, Gai is the Galactic hydrogen column density set equal to 
5.8 x 10 20 cm" 2 (Dickey & Lockman 1990). Using XSPEC 
software (ver. 10.0), we performed a minimum-chi-square 
fitting. The best-fit photon-index is 0.51 yielding x 2 / u = 
210/32 (Figure |l|i). The ratio of the data to the best-fit 
power-law model (Figure Qb) clearly shows evidence for an 
iron K emission line at ~6 keV and a flatter continuum 
above ~8 keV. 

* The pcabackest software was used with the background model 
files of pca_bkgd_faint240_e03v02.mdl and 

pca_bkgd_f aintl7_e03v01 .mdl. 
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Figure 2. The unfolded energy spectrum (detector response was 
approximately removed) of the PCA and HEXTE data (crosses) 
fitted with two power-law model and an emission line. Each model 
component is shown with dotted lines and the total model is the 
solid line. 



Table 1. The RXTE spectral fit 



Figure 1. (a) The background-subtracted spectrum obtained 
with PCA and HEXTE (crosses), fitted with a single power-law 
model convolved with the energy-response matrices (solid lines), 
(b) The ratio between the PCA data and the best-fit power-law 
model. 

Then, we modeled the spectrum with two power-law 
continua (soft and hard) plus a line emission, as expressed by 
f(E) = {K a E- r ° +K h E- Vh +line{K h £ i ))e~ <T P hJVH - Gal . The 
iron- line central energy (Ei) and the intensity (Ki) are free 
parameters, and the line width is assumed to be zero. The 
fit result is shown in Figure ^ and the best-fit parameters 
are summarized in Table 1. Notable properties of the best-fit 
model are very small photon index, ~0, of the hard pawer- 
law continuum, and a large equivalent width, 0.79 keV, of 
the iron line. An absorption edge structure is noticeable at 
7-8 keV in the fit residual (Fig. g), which is most likely 
the iron K-edge. These features are characteristic of the re- 
flected X-rays from an optically thick material as has been 
pointed out by several authors (Iwasawa & Comastri 1998; 
Netzer et al. 1998). X-rays impinging on optically thick mat- 
ter are photo-absorbed as well as Compton-scattered. These 
processes form a very flat continuum around ~ 10 keV, to- 
gether with the K-absorption edge and K emission line of 
iron. We fit the RXTE spectrum together with ASCA data 
in § 4 with models including the Compton reflection. 



3 ASCA DATA 



Parameter Value 



K a (counts s 1 


cm 2 at 1 keV) 


1.83 x icr 2 


r s 




4 4+1-4 


A*h (counts s 


cm - 2 at 1 keV) 


4.45 x 1CT 5 


r h 




04+ 09 


Ki (counts s _1 


cm- 2 )"' 


3.86 x 10" 5 


Ei(keV) 




6.34 ±0.09 


X 




48.0/28 



a > : The physical line width is assumed to be zero. The equiva- 
lent width is calculated to be 0.79 keV. 

3.1 Observation and data reduction 

ASCA observation of NGC 6240 was performed on 27 
March, 1994. The GIS provides an X-ray image with a cir- 
cular field of view of ~ 50' diameter. The SIS was oper- 
ated with 2-CCD mode, covering a much smaller, rectangu- 
lar field of view of ~ 11' x 22'. Iwasawa & Comastri (1998) 
presented the result of the SIS data only. Here, we analyse 
both the SIS and GIS data. In particular, the GIS data are 
useful in order to look for any contamination sources in the 
field around NGC 6240. 

We accepted the GIS and SIS data that were taken when 
the X-Ray Telescope (XRT: Serlemitsos et al. 1995) axis 
was more than 5° above the local horizon, and when the 
geomagnetic cutoff rigidity was larger than 6 GeV/c in order 
to ensure a low and stable background. Additional screening 
condition that the elevation angle from the sunlit earth is 
greater than 25° and 20° was applied to the SIS0 and SIS1 
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data, respectively. As reported by Turner et al. (1997), the 
ASCA data of NGC 6240 show no significant time variation. 



3.2 GIS image 

The GIS image in the full energy band 0.5-10 keV shows a 
point-like source that coincides with NGC 6240 within ~ 1 
arcmin systematic error of the satellite attitude. In addition 
to the emission from NGC 6240, the soft Galactic diffuse 
emission covers the entire GIS field of view. The position 
of NGC 6240, (Z,6) = (20. 73, 27.29), is in the middle of an 
excess-emission structure of Loop I, which is found in the 
ROSAT All Sky Survey (Snowden et al. 1997). The energy 
spectrum of the excess emission is very soft (kT ~ 0.3 keV), 
and it is significant only in the soft energy band below ~ 
2 keV. 

On the other hand, at higher energies, we find an ex- 
tended structure as clearly seen in the image for the 2.4- 
10 keV band shown in Figure in which the background 
consisting of non-X-ray background (NXB) and cosmic X- 
ray background (CXB) has been subtracted. These are es- 
timated from the data obtained in the blank sky observa- 
tions and data taken when the XRT is pointing at the dark 
(night) earth (for detail see Ikebe et al. 1995; Ishisaki et 
al. 1997). The background was subtracted to derive the im- 
age. In Figure ||, notable emission is located north-east of 
NGC 6240, and extends towards north as well as east. Since 
the extension of the emission is much larger than that of the 
soft X-ray emission of NGC 6240 detected with ROSAT, 
we consider it to be an unassociated with NGC 6240. In 
the ROSAT PSPC data, there exist several faint sources 
around NGC 6240 as illustrated in Figure 3. These sources, 
if blurred with the angular response of the GIS, seem to form 
a brightness structure consistent with that observed with the 
GIS. Therefore, the extended X-ray structure is most prob- 
ably due to these faint background sources. Since the field 
of view of RXTE/PCA and HEXTE are both ~1 degree 
FWHM, the emission from the contamination sources will 
contribute to the RXTE spectra significantly. Below we ex- 
amine the ASCA spectrum of NGC 6240 and the properties 
of these contamination sources. 



3.3 Spectrum of NGC 6240 

We construct the GIS energy spectrum of NGC 6240 from 
a circular region of 3 arcmin radius centered on the X-ray 
peak. The spectra from the two GIS sensors, GIS-2 and GIS- 
3, are summed. The NXB + CXB background is subtracted, 
as has been done for the GIS image. The soft excess emission 
is not subtracted, but it is negligible above 2.0 keV. 

The SIS spectrum of NGC 6240 is extracted within a 
circle of 3 arcmin radius from the SIS-0 chip-1 and the SIS- 
1 chip-3. The data taken from the two chips are summed 
together. For the SIS background, we accumulate photons 
from a region where no contamination source is present on 
the same CCD chips with which NGC 6240 was observed. 
Therefore, the soft Galactic diffuse component is also sub- 
tracted as a part of the background. 

The GIS and SIS spectra thus obtained are illustrated 
in Figure [l|. Following Iwasawa & Comastri (1998), we fit 
the ASCA spectrum of NGC 6240 with a two-component 
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Figure 3. The 2.4—10 keV X-ray intensity contour derived with 
the GIS data is shown in logarithmic scale. After the background 
(NXB + CXB) is subtracted, the image was smoothed with a 
Gaussian of a = 1.25'. The outer most solid circle shows the edge 
of the GIS field of view. The position of NGC 6240 is marked 
with *, while • show positions of point like sources detected with 
ROSAT PSPC. 



model that includes a thermal component and an AGN 
component. The thermal component consists of emission 
from two optically-thin thermal plasmas of different tem- 
peratures, where the higher temperature component has an 
excess absorption. The AGN component consists of an ab- 
sorbed power-law continuum and line emission. The model 
can be written as; 



f(E) = {thml(T c ,Z,K c )+e~ 



+ 

x 



line(Ei, Ki) + e 



T ph lv H,AGN 



thml(T h ,Z, K h ) 
KagnE- f } 



(1) 



where the parameters in parentheses are kept free. For the 
thermal plasma emission code, we employ the Mewe-Kaastra 
model (Mewe, Gronenschild, & van den Oord 1985; Mewe, 
Lemen, & van den Oord 1986; Kaastra 1992) modified by 
Liedahl, Osterheld, & Goldstein (1995), which is imple- 
mented in XSPEC as the MEKAL model. The metallici- 
ties, Z, of the two thermal components are assumed to be 
the same. The abundance ratios among different elements 
are fixed to be the fiducial solar values given by Anders & 
Grevesse (1989). 

With this model, we fit the GIS and SIS spectra si- 
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thml(cool) : thml(hot) 

H 1 — I — h- H *f 1 1 H 



power— law 



H 1 1 1 — h- 




0.5 1 2 5 

Energy (keV) 



Table 2. ASCA spectral fit a > 



Parameter 



Value 



K c (cm^ 3 ) 6 ) 

T c (keV) 

K h (cm" 3 ) 6 ) 

T h (keV) 

Z (solar) 

Nn (10 22 cm- 2 ) 

Ki (counts s _1 cm _2 ) c ) 

Ei(keV) 

iV H ,AGN (10 22 cm- 2 ) 
K"agn (counts s — 1 en 

r 

X 2 /u 



0.3V 



3.8 x 10 64 

1+0.11 
-0.05 

3.5 x 10 65 

+0.24 
-0.10 
+ 11.18 
-0.84 
+ 0.14 
-0.16 



1.07 
1.41 
0.92 
2.9 x 10- 



6.46 



+0.03 



0.0 



0.04 
+2.4 



" 2 keV" 1 OlkeV) 2.0 X 10~ 5 

10+ ' 29 
-U.1U_ 43 

59.4/55 



a ) : Throughout this paper, for fitting the SIS spectrum only, an 
additional absorption with equivalent hydrogen column density of 
3.0 X 10 2t) cm -2 was artificially applied, in order to account for a 
systematic overestimate of the absorption column density (Dotani 
et al. 1996; Cappi et al. 1998). 

6 ) : Emission integral. 

c ) : The physical line width is assumed to be zero. The equiva- 
lent width with respect to the power law continuum is 1.24 keV. 



Figure 4. (top) The background-subtracted SIS spectra (crosses) 
fitted with sum of thermal component (dotted lines), an emission 
line (dot-dashed line), and a power-law model (dashed line). All 
the model components are convolved with the energy response 
matrices. The solid line shows the sum of all the model compo- 
nents, (middle) The same as the top panel, but for the GIS data, 
(bottom) The fit residuals for the SIS data (crosses) and the GIS 
data (crosses with an open circle). 

multaneously. In this simultaneous fit, the GIS data below 
2.0 keV are excluded because of a slight uncertainty in the 
energy scale of the GIS arising from the complex xenon M- 
edge structure in this range. A good fit is obtained as shown 
in Figure ^ and Table 2, and the results are essentially the 
same as obtained by Iwasawa & Comastri (1998) using the 
SIS data only. As noted by Iwasawa & Comastri (1998), the 
striking spectral features are a very flat continuum above 
~4 keV with a photon index of ~ and a strong iron K-line 
with an equivalent width of 1.2 keV. 

3.4 Contamination sources' spectrum 

Although the contamination sources near NGC 6240 are not 
identified, we only need the energy spectrum for the purpose 
of the present work. We construct the spectrum from the 
whole field of view of the GIS excluding the 3-arcmin radius 
circle centered on the X-ray peak of NGC 6240, from which 
the estimated background (NXB + CXB) is subtracted. In 
addition, because of an extended outskirts of the XRT point- 
spread function, the spectrum still contains a significant con- 
tribution from NGC 6240. Using a ray-tracing simulation 
software, the contribution from NGC 6240 outside the 3- 
arcmin radius circle is estimated and subtracted. Figure ^ 
shows the GIS spectrum thus obtained. The resultant GIS 



spectrum contains the total photons from the contamination 
sources within the GIS field of view and the Galactic soft 
diffuse emission. 

As shown in Figure H thus obtained GIS spectrum is 
fitted satisfactorily (xV^ — 18.4/25) with the sum of a thin- 
thermal model (MEKAL model) and a power-law model. For 
the thermal component, we obtain the best-fit temperature 
of 0.23±0.04 keV and the surface brightness of 1.0xl0~ 14 
ergs s- 1 cm -2 arcmin -2 (0.5-2 keV), assuming the element 
abundances to be solar. We identify this thermal component 
to be the Galactic diffuse emission, since these values are 
consistent with the ROSAT results (Snowden et al. 1997). 
Thus, the spectrum of the contamination sources can be 
expressed by a power law with the best-fit photon index of 
1.56±0.23. Its flux is 2.1 xlO" 12 ergs s~ x cm" 2 (2-10 keV), 
which is ~ 90% as high as that of NGC 6240. Therefore, the 
contribution from the contamination sources to the RXTE 
spectrum is quite substantial, and will be taken into ac- 
count for the analysis of the RXTE spectrum in the next 
section. It is to be noted that the spectrum of the contam- 
ination sources (added together) is pretty hard, consistent 
with that of AGN. This suggest that these sources, if not 
all, are possibly AGN. 

4 RXTE + ASCA JOINT FIT 

In this section, we model the 0.5-200 keV energy spectrum 
of NGC 6240 obtained with RXTE and ASCA, taking into 
account the contamination sources detected in the GIS field 
of view. Below 3 keV, the X-ray emission is dominated by 
the thermal component that presumably originates from the 
starburst activity in the galaxy, and is well described with 
the two temperature model (Iwasawa & Comastri 1998) 
as shown in § 3.3. The X-ray emission from the AGN in 
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Figure 5. (top) The unfolded energy spectrum of the contamina- 
tion sources obtained with GIS. The spectrum is fitted with sum 
of a thermal and a power-law component (dotted lines). Solid line 
shows the total model, (bottom) The fit residuals. 

NGC 6240 dominates above 4 keV and we perform detailed 
modeling below. 



Figure 6. The result of the simultaneous fit to the ASCA/GIS, 
ASCA/SIS, RXTE/PCA, and RXTE/UEXTE with sum of ther- 
mal component, emission line, and a power-law model, (top) The 
unfolded spectral data (crosses) of ASCA/SIS, RXTE/PCA and 
fiXTE/HEXTE. Individual fitting model components are shown 
in dotted lines, while the total is given with the solid line. The 
GIS data are omitted for the display, (middle) The fit residuals 
for the SIS data (crosses) and for the GIS data (crosses with open 
circle), (bottom) The fit residuals for the PCA and HEXTE data. 



4.1 Power law model 

For modeling the energy spectrum of the AGN, we begin 
with the simplest model, power-law + line, given by equa- 
tion |l| in § 3.3, The ASCA GIS and SIS spectra of NGC 6240, 
the PCA spectrum, and the HEXTE spectrum were fitted 
simultaneously. Since the X-ray intensity of the AGN may 
have varied between the ASCA and RXTE observations per- 
formed at different periods, the normalization factor of the 
power-law component, T^agn, is left free for the ASCA and 
RXTE spectra, respectively. Since the iron K-line is pre- 
sumably generated by reprocessing of the X-rays from the 
AGN on an optically thick matter, the line intensity will 
be proportional to the AGN luminosity on the time aver- 
age. Therefore, we assumed that the ratio of the flux of the 
iron line to that of the power-law continuum, Ki/Kagn, was 
the same in the two separate observations with ASCA and 
RXTE, respectively. On the other hand, the normalization 
factor of the thermal components, K c and Kh , are assumed 
to be common to the RXTE and ASCA spectra, since the 
thermal component which is extended should not vary with 
time. All other parameters are also tied between the two 
observations. In order to take into account the hard con- 
tribution source, a power-law model with the photon index 
1.56 as derived in the previous section is included in the 
PCA and HEXTE spectra with a free normalization factor. 



We assume that the power law spectrum extends beyond 
10 keV with the same photon index. 



As shown in Figure |fj and Table 4T , the fit is not accept- 
able. As expected, the obtained parameters of the thermal 
component are essentially the same as those given in Sec- 
tion 3.3, since the RXTE spectrum has little influence below 
4 keV. However, the best-fit power-law continuum of AGN is 
significantly flatter than that given by Iwasawa & Comastri 
(1998). Furthermore, the residuals require an even harder 
continuum in 8-15 keV band, and also suggests a spectral 
steepening above ~ 30 keV. These features together with 
a large equivalent width of Fe K-line are consistent with 
a Compton reflection spectrum. We therefore incorporate 
Compton reflection in the next subsection. 



4.2 Reflection model 

We employ the Compton reflection model developed by 
Magdziarz & Zdziarski (1995) implemented in XSPEC as 
the PEXRAV model. Using the PEXRAV model, we can cal- 
culate the reflected X-ray spectrum when X-rays of a power- 
law spectrum illuminate an optically thick layer of material 
that is predominately neutral except hydrogen and helium. 
Here, we assume that the direct X-rays from the AGN are 
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Table 3. Power law model 



Parameter 


Value 


K,(ASCA)(counts s" 1 cm- 2 )"' 


2.8 X icr 5 


£,(keV) 


6.44 


^H.AGN (10 22 cm- 2 ) 


1.1 


r 


-0.15 


^agn(ASCA) 6 ) 


1.9 x 10~ 5 


K AGN (RXTE)/K AGN (ASCA) 


1.16 


X 2 /" 


131.5/87 



a ) : The physical line width is assumed to be zero. The equiva- 
lent width with respect to the power-law continuum is 1.13 kcV. 

b ) : In the unit of counts s" 1 cm" 2 keV" 1 CSlkeV. 

totally blocked by a thick absorber on the line of sight. The 
parameters that describe the model are: V, photon index of 
the incident power law spectrum with an exponential cutoff 
at E c ; R = Q/2tt, where Q is the solid angle subtended by 
the optically thick material; /i = cos8, where 8 is the an- 
gle between the line of sight and the normal vector of the 
optically thick layer; A-p e , the iron abundance of the opti- 
cally thick material. Element abundances besides iron are 
assumed to be the solar values. 

In the following fitting, E c is fixed at 200 keV (practi- 
cally a single power law in the observed range), /j, is fixed 
at 0.45 so that the model spectrum is closest to that aver- 
aged over all viewing angles as described in Magdziarz & 
Zdziarski (1995). Combining with the thermal component 
and the line emission, the overall fitting model is given by 

{thml(T c , Z, K c ) + e-^^H thm i(r h , Z , K h ) 



f(E) 



+ line(Ei,Ki) + refl(r,K A GN)}e- 



[T ph- iV H,Ga,l 



(2) 



where free parameters are shown in parentheses. The nor- 
malization factors of the reflection component, -Kagn, for 
the ASCA and RXTE spectra are left free from each other 
for possible time variation of the AGN luminosity. The 
flux ratio between the Fe K-line and the AGN continuum 
Ki I A"agn is assumed to be constant. All other parameters 
are also assumed to be the same between the two observa- 
tions. Another power-law with a photon index of 1.56 that 
represents the contamination sources is added to the model 
for the PCA and HEXTE spectrum. 

As shown in Figure Q and Table 4, this model can ac- 
count for the 0.5-200 keV wide-band spectrum satisfactorily. 
The flux of the contamination sources in the RXTE spec- 
trum determined from the fit is in good agreement with that 
obtained with the GIS within ~ 10%. The two normaliza- 
tion factors of the AGN component obtained from the ASCA 

and RXTE spectra respectively happened to be very close 
^RXT E / jr ASCA =0 92 _ 1 26 ^ ag were the Fe K _ line intensi . 

ties nearly equal in these two observations. The equivalent 
width of the iron line with respect to the reflection contin- 
uum is 0.8 keV, which is a typical value for the fluorescence 
line associated with the Compton-reflection. 

However, the best-fit photon index P of the incident 
AGN spectrum is 1.26±0.13, which is unusually small com- 
pared to the typical values for AGNs (e.g. Mushotzky, Done, 
& Pounds 1993). In this fit, the viewing angle /i = cos9 was 
fixed at 0.45. Even if we allow /i to vary (Table 4), photon 




20 50 
Energy (keV) 



200 



Figure 7. Unfolded spectrum of the RXTE data (crosses) fitted 
with sum of the thermal component, reflection component, and a 
power-law for the contamination source. Each model component 
is shown with dotted lines, while the total is illustrated with the 
solid line. Although the ASCA/GIS and SIS data are not shown 
here, they are also fitted jointly. 



index larger than 1.54 is not acceptable. The intrinsic lumi- 
nosity of the AGN is estimated to be larger than 4.3 x 10 43 
ergs/s, the value corresponding to the maximum solid angle, 
i.e. R = n/2-K, of 1.0. 



4.3 Reflection + Transmission model 

Although, the above reflection-only model gives a satisfac- 
tory fit to the RXTE + ASCA data, the small photon index 
derived in § 4.2 still remains to be problem. As shown below, 
this is resolved by adding an absorbed AGN continuum that 
is transmitted through a thick absorber on the line of sight. 
(Here we do not consider the extent of the absorber, hence 
we assume no scattering of X-rays into the line of sight.) 
The transmitted AGN component is given by 

-E/E c 



(3) 



where cttii is the Thomson scattering cross-section, and 
A'h.agn is the hydrogen-column density along line of sight. 
The attenuation by the Thomson scattering becomes impor- 
tant above ~ 10 keV as compared to photoelectric absorp- 
tion. 

Adding the transmitted AGN component given by eq. ^ 
to eq. g, we construct the reflected- and transmitted-AGN + 
thermal model as; 



f{E) 



{thml{T c ,Z,K c ) 



7 ^ Nh thml(T h ,Z,K h 
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Table 4. reflection model 



Parameter 



Value (fit l) a ) Value (fit 2) a > 



r 


1 26+ - 13 


0.98±g : » 


H = cosO 


0.45 (fix) 


«o.08lg;S 


Fe abundance (solar) 


n 47+O.22 
u '^' -0.19 


n yo+O.SS 
U -' Z -0.32 


^agn(ASCA) c > 


2.7xl0~ 3 


4.6xl0~ 3 


Ki{ASOA) (counts s" 1 cm- 2 ) d ) 


2.5xl0" 5 


2.4xl0" 5 


E l (ASCA)(kcV) 




6.44±0.04 


K AGN (RXTE)/K AGN (ASCA) 


i-08±8:1S 


1 05+ 016 


X 2 /» 


97.2/87 


96.4/86 



a > '. fi is fixed or left free in the fit 1 and 2, respectively. 

: The lower and upper limits correspond to the limitation of 
the model code. 

c ' : In the unit of counts s~ 1 cm -2 keV -1 @lkeV. 

d ) : The physical line width is assumed to be zero. Equivalent 
width with respect to the reflection continuum is 0.88 keV and 
0.84 keV for the fit 1 and 2, respectively. 



> o 
M 



si 
IX 



X o 
I 




10 



20 50 
Energy (keV) 



100 



200 



+ line(Ki,Ei) + refl(T,R,E c ,K AGN ) 

+ e -N R , AGN ^ h +a Tb ) KAG ^ E -V e -B/B c y 

X g _ °'ph JV H,Gal 



(4) 



The photon-index F, the cutoff energy, Ec (fixed to 
200 keV), and the normalization, T^agn, are the same for 
both the reflection and transmitted components. The rela- 
tive intensities between the transmitted component and the 
reflection component is determined by the solid angle, R, 
which is left free. The metallicity of the elements lighter 
than Fe and the viewing angle /i = cos9 are fixed at 1.0 so- 
lar and 0.45, respectively. As in the previous fits, including 
the contamination-source component in the RXTE model, 
we fit the four spectra jointly with this model. 

The reflected- and transmitted-AGN + thermal model 
gives a good fit (Figure ^, Table 5). The derived photon 
index is F= 1.59 io' 4 <3' wmcn ^ s consistent with the canonical 
value of 1.9 (Pounds et al. 1990). The best-fit solid angle of 
the reflector, R = Q/2ir, is 0.51, and the absorption column 
density, A?h,agn, is found to be 1.7 x 10 24 cm -2 . Since R 
was poorly constrained, the intrinsic luminosity of the AGN 
has a large error range. If we assume that R does not exceed 
1.0, I/x is obtained to be l.llg'5 x 10 44 ergs/s in the range 
2-10 keV. When we leave the viewing angle fi free, the 
fitting parameters as well as the intrin sic l uminosity remain 
essentially the same, if R < 1.0 (Table 4.3). 

So far, we assumed that the matter that obscures the 
AGN is located only on the line of sight. If the absorbing 
matter covers a significant solid angle viewed from the AGN 
such as in the case of a torus, some of the incoming photons 
will be scattered into our line of sight. Then the model given 
by eq. |§] would overestimate the AGN luminosity. Matt, 
Pompilio, & Franca (1999) performed a Monte Carlo cal- 
culation of the X-ray transmission through spherically dis- 
tributed matter for various column densities. According to 
their result, the true AGN luminosity would be smaller by 
a factor of ~ 2, if the absorbing matter is distributed in 
a spherical geometry. Consequently, for the same flux from 
the cold reflector, the true value of R should be larger by a 



Figure 8. The fit result of the reflection-and transmitted- AGN 
model. Crosses are the unfolded spectral data. Each model com- 
ponent is shown with dotted lines, while the total is given with 
solid line. Although the ASCA/GIS and SIS data were also fitted 
jointly, only the RXTE data are shown here. 



factor of 2 than that in Table 5. The column density of the 
absorbing matter is also subject to a slight overestimation. 

The case that a heavy absorber exists only on the line of 
sight is rather unlikely. On the other hand, since the fitting 
result shows a relatively small absorption for the reflection 
component, the absorbing matter is probably not covering 
the entire sphere. There two cases, a cloud on the line of 
sight are spherically distributed matter, are considered to 
represent two extremes. 

In conclusion, the reflected- and transmitted-AGN 
model gives the following AGN parameters: F = 1.33 — 2.02, 
L x (2-10keV) = 5 x 10 43 - 6 x 10 44 ergs/s, R = 0.1-1, and 
A^h.agn = 1.0 - 2.7xl0 24 cm" 2 , where the range of param- 
eters represents not only the statistical errors but also the 
uncertainties due to unknown inclination angle of the reflec- 
tor, and the limits for the two extreme geometries of the 
absorbing/scattering material. 



5 DISCUSSION 

The X-ray spectrum of NGC 6240 is satisfactorily explained 
with a model consisting of a thermal component and an 
AGN component. For the AGN component, and equally 
good fit is obtained either by the reflection-only model or by 
the reflected- and transmitted-AGN model. Hence, we can- 
not conclude the detection of the direct AGN component 
from the present data. However, without the direct (trans- 
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Table 5. Reflection and transmitted AGN model 



Parameter 


Value (fit 1) 


Value (fit 2) a > 


r 


i.oj_ 26 


, 71 +0.3l 
1 - ' -0.37 


H = cos9 


0.45 (fix) 


6) gl +0.04 


Fe abundance (solar) 




fl 77+1-23 


^agn(ASCA) c > 


9.1xl0~ 3 


1.3xl0" 2 


A^agn (xl0 24 cm- 2 ) 


1 -' -0.5 


i 7+O.7 
1 -'-0.5 


R 


(1 51+°" 
U - Oi _0.41 


fl q^+0 66 a) 
U - J *-0.21 


ii";(counts s~ 1 cm" 2 ) 1 *' 


2.4xl0" 5 


2.3xl0" 5 


Ei(keV) 


6 44+ 04 
°- 44 -0.05 


6 44+ 04 
Oi44 -0.05 


K AGN (RXTE)/K AGN (ASCA) 


1 00+"' 18 


1 03+ 018 


X 2 /v 


86.2/85 


84.8/84 



a ' : In the fit, fj, is left free. Maximum value of R is assumed to 
be 1. 

" : The lower and upper limits correspond to the limitation of 
the model code. 

c ) : In the unit of counts s _1 cm" 2 keV -1 CSlkeV. 

d ) : The physical line width is assumed to be zero. The equiv- 
alent width with respect to the reflection continuum is 0.96 keV 
and 0.95 keV for the fit 1 and 2, respectively. 



mitted) component the derived photon index is noticeably 
smaller than the typical values for AGN. We consider it more 
probable that the direct AGN component is present at high 
energies. 

According to the unified model for Seyfert galaxies (e.g. 
Antonucci 1993), one can interpret the result in terms of a 
tilted molecular torus in which the near-side of the torus 
acts as an absorber and the far-side acts as a reflector. If 
the torus is axisymmetric, both the reflecting part and the 
absorbing part would have a similar column density. With 
the reflected- and transmitted-AGN model, the column den- 
sity of the obscuring matter is estimated to be ~ 2 x 10 24 
cm~ 2 , corresponding to Thomson optical depth of ~ 1. This 
is practically enough for a reflector. Thus, the present result 
is consistent with the interpretation that NGC 6240 has a 
similar geometry to Seyfert 2 galaxies. The obtained column 
density gives an IR extinction of Ak ~ 100, which is consis- 
tent with the fact that the previous IR observations did not 
find the AGN itself. Narrow line region that has not been 
seen in the optical observations would be obscured as well. 

The intrinsic X-ray luminosity of NGC 6240 is esti- 
mated in § 4.3 to be in the range 5 x 10 43 — 6 x 10 44 ergs/s 
in the range 2-10 keV. For the reflection-only case (§ 4.2), 
the luminosity is estimated to be larger than 4 x 10 43 ergs/s, 
while the upper bound of luminosity is not determined. How- 
ever, it would be plausible that the solid angle factor R is not 
much less than several %, which gives essentially the same 
luminosity range as the above. This luminosity is among 
those of the most luminous Seyfert nuclei, and even compa- 
rable to those of quasars. The bolometric luminosity of the 
AGN may well exceed 10 45 erg/s. 

Concerning the power source of the huge IR lumi- 
nosity, a measure of the contribution of an AGN to the 
IR emission is a ratio of X-ray luminosity to IR luminos- 
ity, Lx/Lxr. Based on our results given above and the 
IR flux calculated from the IRAS result with the formula, 



Fir = flux(25^m) X (^25 M m) + flux(60^m)x(i'60Mm), gives 
L x (2 - 10 keV)/L m = 0.01-0.1 for NGC 6240. This agrees 
with those of other Seyfert nuclei given by Ward et al. 
(1988). 

Vignati et al. (1999), using the BeppoSAX data of 
NGC 6240, claim that the direct AGN component was pos- 
itively detected, though they assumed the photon index of 
the AGN to be 1.8. Other than that, their results and inter- 
pretation are essentially in agreement with ours. 

The luminosity of the AGN can be used to infer the 
mass of the central black hole in NGC 6240. If the AGN 
accretes about 1% of the Eddington luminosity, the total 
luminosity of 10 45 erg/s would lead to a mass of Magn — 
10 9 Mq, which is consistent with the value expected from the 
galaxy-mass to black-hole-mass relation (Lauer et al. 1997), 
with the estimated galaxy mass of 10 11_12 M Q after having 
completed its merging epoch (Shier & Fisher 1997). 



6 CONCLUSION 

The 0.5-200 keV wide band energy spectrum of NGC 6240 
obtained with RXTE and ASCA is accounted for in terms 
of a soft thermal component and a hard AGN component. 
The soft component is presumably due to star burst activ- 
ity. The AGN component consists of a Compton reflection 
component accompanied by an intense Fe-K emission line 
and probably a transmitted component (a direct component 
penetrating through a thick absorber) . The detection of the 
transmitted component is not conclusive from the fitting. 
However, without a transmitted component, the photon in- 
dex is unusually small. Assuming that the solid angle factor 
R does not exceed 1, we estimated the intrinsic X-ray lumi- 
nosity of the AGN in the range 2-10 keV to be in the range 
4 x 10 43 — 6 x 10 44 ergs/s, which yields the ratio of the X-ray 
luminosity (2-10 keV) to the IR luminosity of 0.01-0.1. The 
column density that obscures the central AGN is estimated 
to be larger than 1.0 x 10 24 cm -2 . These results show that 
NGC 6240 is among the most luminous Seyfert 2 galaxies. 
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